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Abstract. With Monte Carlo methods we study the dynamic relaxation of a vortex 
state at the Kosterlitz-Thouless phase transition of the two-dimensional XY model. 
A local pseudo-magnetization is introduced to characterize the symmetric structure 
of the dynamic systems. The dynamic scaling behavior of the pseudo-magnetization 
and Binder cumulant is carefully analyzed, and the critical exponents are determined. 
To illustrate the dynamic effect of the topological defect, similar analysis for the the 
dynamic relaxation with a spin-wave initial state is also performed for comparison. 
We demonstrate that a limited amount of quenched disorder in the core of the vortex 
state may alter the dynamic universality class. Further, theoretical calculations based 
on the long-wave approximation are presented. 



PACS numbers: 64.60.Ht, 68.35.Rh, 05.10.Ln 



Dynamic relaxation of topological defect at Kosterlitz-Thouless phase transition 



2 



1. Introduction 

In recent years many activities have been devoted to the study of dynamic processes 
far from equilibrium. Compared with spin-glass and structural-glass dynamics, critical 
dynamics at standard second-order or Kosterlitz-Thouless (KT) phase transitions is 
relatively simple. One is able to systematically explore the universal dynamic scaling 
behavior far from equilibrium, up to the macroscopic short-time regime [TJ [2j El 0]. 
Although the spatial correlation length is still short in the beginning of the time 
evolution, the dynamic scaling form is induced by the divergent correlating time. Based 
on the short-time dynamic scaling form, new methods for the determination of both 
dynamic and static critical exponents as well as the critical temperature have been 
developed [3j HJ E] . Since the measurements are carried out in the short-time regime, 
one does not suffer from critical slowing down. Such a dynamic approach is especially 
useful, when one is interested in tackling both the statics and dynamics of the critical 
systems. Recent progress in the short-time critical dynamics includes, for example, 
theoretical calculations and numerical simulations of the XY models and Josephson 
junction arrays [HI El EJ [9j HO], magnets with quenched disorder (TTJ [12} [13j [HI [15], 
ageing phenomena [H3, [T71 UHl HH ED], domain- wall dynamics [211 E21 [231 EI], weak 
first-order phase transitions [251 EH1 [HI EI] , and various applications and developments 

[2H1 EHl EQl EH E21 [331 EH [35]. 

In the understanding of the critical dynamics far from equilibrium, the dependence 
of the dynamic scaling behavior on the macroscopic initial condition [HE! [36] is essential. 
Up to now, the dynamic relaxation with disordered and ordered initial states has been 
systematically investigated [U El El [36]. Physically, the disordered and ordered initial 
states correspond to the states at very high and zero temperatures respectively. The 
ordered initial state can be also considered as a state under a strong external magnetic 
field. Recently, the dynamic relaxation of a domain wall is also concerned [2H 122] . 
and it shares certain common features with those around free and disordered surfaces 
|21j . The three relaxation processes above are typical but relatively simple. By these 
successes, we are encouraged to tackle more general relaxation processes. On the other 
hand, for strongly disordered systems such as spin glasses, not only the metastable states 
but also the ground states are not simply ordered and homogeneous, and the standard 
magnetization does not characterize the relaxation dynamics starting from a zero or low 
temperature. Methodologically and technically one needs to develop new concepts. 

From theoretical view points, metastable states and low-energy excited states may 
play important roles in phase transitions. One example is the spin-glass transition, and 
another is the KT transition in the two-dimensional (2D) XY model. In the latter case, 
low-energy excited states such as spin waves destroy the magnetic order, and metastable 
state such as vortices and vortex pairs dominate the phase transition. Since vortices and 
vortex pairs are topological defects, the KT transition is sometimes called a topological 
phase transition. The effect of topological defects is often a focus issue in the study of 
phase transitions, also in non-equilibrium dynamics [37] . In fact, one has not achieved 
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too much understanding for the dynamic relaxation of metastable states and low-energy 
excited states, although its importance is often qualitatively addressed. 

Recently many activities have been devoted to the study of vortex states and vortex 
structures in nanowires and nanomagnets. Attention is also drawn to the single vortex 
state in experiments and numerical simulations [38l [39l EH EJ H2] . For example, a single 
vortex structure is explicitly observed in the domain wall of a ferromagnetic nanowire 
[40J, and a second-order transition between a domain wall and a vortex state is detected 
[42] . In Ref. [41] . the abnormal diffusion of a single vortex in the 2D XY model is also 
numerically simulated. Therefore, the vortex state is a physically important object. In 
this paper, we numerically and theoretically study the relaxation dynamics of a single 
vortex state at the KT phase transition of the 2D XY model, in comparison with that of 
a spin-wave state. The dynamic effect of the topological defect will be emphasized. We 
demonstrate that it is essential to introduce a local pseudo- magnetization to characterize 
the symmetric structure of the dynamic systems, which is defined as the projection of 
the spin configuration onto the metastable state or low-energy excited state. In Sec. 
2 the model and scaling analysis are described, and in Sec. 3 the numerical results 
are presented. In Sec. 4 theoretical calculations based on the linearized long-wave 
approximation are presented. Sec. 5 includes the conclusions. 

2. Model and scaling analysis 

2.1. Model 

The 2D XY model is the simplest model for magnetic materials, exhibiting a KT phase 
transition. The Hamiltonian is written as 

K1 <ij> 

where S% = {Si >x , S^y) is a planar unit vector at site % of a square lattice, the sum is 
over the nearest neighbors, and T is the temperature. In our notation, we simply take 
K = 1/T. In the literatures, the transition temperature T KT is reported to be between 
0.89 and 0.90 for the 2D XY model [43j HU H5] . Below Tkt, the system remains critical. 
To study the pure relaxation dynamics of model A, we adopt the "heat-bath" algorithm 
of an one-spin flip, in which a trial move is accepted with probability l/[l+exp(A£'/T)], 
where AE is the energy change associated with the move. 

To study the dynamic relaxation of a single vortex state, we first construct a perfect 
vortex as the initial state. We put the vortex on a square lattice L x L, and choose 
the center of the vortex as the origin of the polar coordinates. To reduce the finite- 
size effect, we apply antisymmetric boundary conditions with respect to the origin, i.e., 
S? = —£>-r, which mimic the vortex line on the boundary. In Fig. Ufa), the dynamic 
evolution of the spin configuration is illustrated. Due to the vortex initial state, the 
dynamic system is inhomogeneous for different radius r. After the statistical average, 
the orientations of the spins form circular vortex lines. Except for the topological defect 
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at the center, the vortex state is locally ordered, although globally it is still different 
from an ordered state. Naively, we may measure the magnetization along the circular 
vortex lines as a function of r and t. From the symmetry of the initial state, however, 
both the x and y components of the magnetization are zero. 

In order to characterize the dynamic relaxation of the vortex state, therefore, we 
should introduce a local pseudo-magnetization and its second moment by the projection 
to the metastable state, 




N r 



<p=l 



rip,p 




S ri p{t) ■ s. 

<p=l 



r<p,n 



k = 1,2; 



fc = l,2, 



(2) 



(3) 



where (■ • ■) represents the statistical average, and N r w 2nr is the number of the lattice 
sites at radius r. In the case now, the metastable state is just the perfect vortex, 



the same as the initial state, S, 



rip,p 



S rtp (0) 

that of the normal direction. 



e v , while S r(Ptn 



e r denotes the 
In other words, 



unit vector perpendicular to S r(PiP , i.e 

r) is the k-th moment of the tangent direction, and M^ k \t,r) is that of the 
normal direction. From the symmetry of the vortex initial state, the normal component 
of the pseudo-magnetization is zero, i.e. M n (t,r) = 0. Denoting M(t,r) = M p (t,r) and 
M^ 2 '(t,r) = Mj®(t,r) + Mj®(t,r), we then define a time-dependent Binder cumulant, 

M {2 \t,r 



U(t,r) 



- 1. 



(4) 



M(t,r) 2 

The Binder cumulant U(t, r) describes the fluctuation of the pseudo-magnetization. 

To investigate the dynamic effect of quenched disorder, we may randomly fix the 
orientations of the four spins at the center of the vortex during the dynamic evolution. 
Such a dynamic system is topologically similar to the dynamic relaxation around a 
disordered surface [21], [31]. Experimentally, a vortex state may be prepared by applying 
an electric current perpendicular to the lattice plane. 

For comparison, we now consider the dynamic relaxation starting from a low-energy 
excited state, i.e., a spin- wave state. The spin wave is linearly polarized in the x 
direction, while ordered in the y direction. The simulation is performed with a rectangle 
lattice L x xL y . In Fig.[T](b), the dynamic evolution of the spin configuration is displayed. 
The spin- wave state is locally ordered, but globally not. For such a dynamic process, 
we may naively measure the line magnetization along the x direction as a function of t. 
Obviously, the magnetization is zero due to the periodicity of the spin-wave state. 

In order to characterize the dynamic relaxation starting from the spin-wave state, 
we may also define a pseudo-magnetization and its second moment as functions of t, 



Mf (t) 



Tk 



xy.p 



k = 1,2: 



(5) 
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xy,n 



k = 1,2, (6) 



where = •S'xy(O), denotes the unit vector perpendicular to S X y,p, and (■ ■ ■) 

represents both the statistical average and the average in the y direction. As shown 
in Fig. QJb), S xy (0) periodically changes its orientation along the x direction, but 
independent of y. From the periodicity of the spin- wave state, M n (t) = 0. Denoting 
M(t) = M p (t) and M^{t) = Mj®(t) + M®(t), we then define a time- dependent Binder 
cumulant, 

2.2. Scaling analysis 

In the critical regime, there are three spatial length scales in the dynamic relaxation 
of the vortex state, i.e., the nonequilibrium spatial correlation the radius r of the 
circular vortex line, and the lattice size L. Therefore general scaling arguments lead to 
the scaling forms of the magnetization and its second moment, 

M (i) (i, r, L) = at)~ kvl2 M {k \r/m,m/L), k = 1, 2, (8) 

where 77 is the static exponent. In the equation, the factor £ > {t)~ kr) l 2 indicates the scaling 
dimension of M^, and the scaling function M^(r/£(i), £(i)/L) represents the scale 
invariance of the dynamic system. In general, we expect that the scaling forms in Eq. (jSJ) 
holds already in the macroscopic short-time regime, after a microscopic time scale t m i C 
[U [3] . Our numerical simulations show, however, that only the pseudo-magnetization 
and its second moment obey such scaling forms. 

In the short-time regime, i.e., the regime with <C L, the pseudo-magnetization 
is independent of L. Then the scaling form is simplified to 

M(t,r)=at)-" /2 M(r/at)). (9) 

Since the spatially correlating terms in the susceptibility M^ 2 ^(t,r) — M(t,r) 2 can be 
neglected for £(£) r, one may deduce U(t,r) ~ for large r, with d = 2. Here 

N r oc r is the number of the lattice sites at radius r, similar to the lattice size L in 
Refs. [31 |2U [22]. Further, U(t,r) is independent of the scaling variable £,{t)/L in the 
short-time regime. Together with Eqs. (IE]) and (Q, one may write down the scaling form 

U(t,r) = m/r] d - 1 U(r/m)- (10) 

The Binder cumulant is interesting, for the static exponent 77 is not involved. 

For the critical dynamics of a continuous phase transition, £(t) usually grows by a 
power law f (t) ~ t x l\ and z is the so-called dynamic exponent. In shorter times, there 
may be corrections to scaling, typically in a power-law form 

m-t 1/z a+c/t b ). (11) 
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For magnetic systems with a KT phase transition, e.g., the 2D XY model, the correction 
to scaling is weak in the dynamic relaxation with an ordered initial state, the dynamic 
exponent is theoretically expected to be z = 2 and numerical simulations indicate b ~ 1 
[HI [19]. Due to the dynamic effect of the vortex-pair creation and annihilation, however, 
the correction to scaling becomes strong in the dynamic relaxation with a disordered 
initial state, and both theoretical and numerical calculations lead to a logarithmic form 

[an m ei mi us] 

m-WVnt + c)) 1 /*. (12) 

Theoretically, Eq. ( 112]) is equivalent to Eq. (II ip in the limit b — > 0. Numerically detecting 
a logarithmic correction to scaling is rather notorious, for it is negligible only in the limit 
t — > oo [6J 113 [19]. Recently, a logarithmic correction to scaling is detected inside the 
domain interface, which is also attributed to the vortex-pair creation and annihilation 
[22]. 

In this paper, we will show that there exists a core of the vortex state, and there 
emerges a strong logarithmic correction to scaling in the core, again due to the vortex- 
pair creation and annihilation. The core of the vortex grows with time, and also 
roughens. Such a phenomenon is rather similar to the propagation and roughening 
of the domain interface [211 [221 [23]. Outside the core of the vortex and also for the 
spin-wave initial state, the dynamic system is locally ordered. Therefore, the dynamic 
relaxation is similar to that with the ordered initial state, and the correction to scaling is 
in a weak power-law form. Further, we demonstrate that a limited amount of quenched 
disorder in the core of the vortex may change the dynamic universality class, i.e., the 
scaling functions M^(s) in Eq. ([HJ) and the corresponding critical exponents. 

For the dynamic relaxation with the spin-wave initial state, there are two spatial 
length scales in the dynamic system, i.e., the nonequilibrium spatial correlation £(£), and 
the lattice size L x , if the lattice size L y is sufficiently large. Therefore scaling arguments 
lead to the scaling form of the pseudo-magnetization and its second moment, 

M( k \t,L x )=Z(ty kv/2 M( k) m/L x ), k = l,2. (13) 

For the pseudo-magnetization, the scaling function M(^(t)/ L x ) is independent of L x in 
the thermodynamic limit L x — > oo. The scaling form is simplified to 

M(t) ~ i{t)^ /2 . (14) 

The scaling behavior of the Binder cumulant is, 

U(t) ~ m/L x ] d - 1 - (15) 



3. Monte Carlo Simulation 



For the dynamic relaxation with the vortex initial state, our main results are obtained 
with L = 256 at Tkt = 0.89, and the maximum updating time is tu = 10240 Monte 
Carlo time steps. A Monte Carlo time step is defined by a sweep over all spins on the 
lattice. Additional simulations with L = 128 and L = 512 are performed, to investigate 
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the finite-size scaling behavior and finite-size effect. Total samples for average are about 
20000. For the dynamic relaxation with the spin-wave initial state, the main results are 
obtained with L x x L y = 256 x 256 and L x x L y = 512 x 256 at T KT = 0.89, and the 
maximum updating time is tu = 10240 Monte Carlo time steps. Additional simulations 
with L x xL y = 128 x 256 and L x xL y = 1024 x 256 are also performed. Total samples for 
average are about 10000. Statistical errors are estimated by dividing the total samples 
into two or three subgroups. If the fluctuation in the time direction is comparable with or 
larger than the statistical error, it will be taken into account. Theoretically, the scaling 
forms described in the previous section hold in the macroscopic short-time regime, after 
a microscopic time scale t m i C . t m i C is not universal, and relies on microscopic details of 
the dynamic systems. In our simulations, t mic is about 100 Monte Carlo time steps. 



3.1. Vortex state 

For the dynamic relaxation with the vortex initial state, the time evolution of the local 
pseudo- magnetization M(t,r) = M^'(t,r) defined in Eq. fl2]) is displayed in Fig. |2]^a). 
Monte Carlo simulations have been performed with the lattice sizes L = 256 and 512, 
and the finite-size effect is negligibly small. Denoting s = r/£(t), we observe that for a 
sufficiently large s, e.g., r = 127.5 and t < tu = 10240, M{t,r) approaches the power- 
law decay at bulk, i.e., M(t,r) ~ t~ v ^ 2z , while for a sufficiently small s, e.g., r = 0.707 
and t > 100, M(t, r) also exhibits a power-law behavior, but decays much faster than 
at bulk. In other words, the scaling function M(s) in Eq. ([9]) is characterized by 

~, , r const s — > oo ,.„ s 

s ^o ' (16) 

Here rjo is the vortex exponent. We call the regime with small s the core of the vortex 
state. The spatial length scale of the core increases with time, and it is just proportional 
to the nonequilibrium correlation length £(£). Outside the core of the vortex, i.e. the 
regime with large s, the dynamic relaxation of the pseudo-magnetization is governed by 
the bulk exponent 77, while in the core of the vortex, it is controlled by both the bulk 
exponent n and the vortex exponent r)o- Outside the core of the vortex, the dynamic 
relaxation of the pseudo-magnetization is the same as that with an ordered initial state, 
and the correction to scaling is weak. Assuming £(£) ~ t 1 ^, one deduces M(t, r) ~ t~ v ^ 2z 
for s — > 00. In Fig. |2^a), the exponent t]/2z = 0.0592(4) measured from the slope of 
the curve of r = 127.5 is in agreement with that for the ordered initial state (6j [22] . 
and it is also consistent with rj = 0.234(2) and z = 2 reported in the literatures [6]. 
A power-law correction to scaling slightly refines the result by about one percent, and 
yields t]/2z = 0.0586. In fact, the exponent 77 extracted in almost all the numerical 
simulations is somewhat smaller than the prediction of the KT theory, rj = 1/2. This 
may imply that the KT theory is not exact for the lattice XY model. In the core 
of the vortex, one may similarly derive M(t, r) ~ 

t -(r,+ m )/2z for s _^ q_ Then 

one 

measures (t] + 7] )/2z = 0.498(4) from the slope of the curve of r = 0.707, and calculates 
?7o/2 = 0.878(8) by taking z = 2 as input. 
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For the domain interface of the 2D XY model, the interface exponent 770/2 is 
reported to be 0.997(7), very close to 1 [22J. It indicates that M(t,r) is an analytic 
function of r. Topologically, a vortex state could be transformed from either a domain 
interface or a free surface. Physically, the former seems closer to a vortex state. 
Is the vortex exponent r/o/2 really different from 1? According to the scaling form 
in Eq. ffTB"]) . the r-dependence and t-dependence of the magnetization M(t,r) should 
yield a same exponent 770/2. We observe that at a fixed time t, the r-dependence 
of M(t,r) is very close to linear in the small r regime. According to Eq. ( fT6l) . this 
indicates 7/0/2 rs 1, contracting to r/ /2 = 0.878(8) from the power-law fit in Fig. [2(a)- 
Our thought is that there exists a strong logarithmic correction in the growth law of 
£(t), described by Eq. f|T2l . Such a fitting to the curve of r = 0.707 in Fig. [2(a) 
yields (77 + r] )/2z = 0.550(9). With rj/2z = 0.0586 as input, one then calculates the 
vortex exponent 770/2 = 0.983(18), also close to 1. To further confirm this result, 
let us fully examine the scaling form in Eqs. (Q and (|T6l) . In Fig. 12(b), the scaling 
function M(r/£(t)) = M(t, r)^{t) r] / 2 is plotted, and data collapse is clearly observed for 
different time t. In the small-s regime, the logarithmic correction to scaling described 
by Eq. ( jT2l) has been taken into account, and the constant c = 5.45 is taken from the 
fitting in Fig. [2(a). Obviously, M(s) exhibits a power-law behavior in the small- s regime, 
supporting the linear r-dependence of M(t, r). From the slope of the curve of t = 10240, 
one measures ?7o/2 = 0.99(2), rather close to 1, and consistent with 770/2 = 0.983(18) 
obtained from Fig. [2(a). For large s, M(s) approach a constant. Finally, we mention 
that in the data collapse in Fig. [2(b), data points for all r including r = 0.707 are 
plotted, although those for t < t mic ~ 100 or 200 are truncated. In other words, the 
dynamic scaling behavior holds even in the very center of the vortex state. 

The logarithmic correction to scaling at the KT phase transition is believed to be 
induced by the vortex-pair creation and annihilation. Therefore we measure the time 
evolution of the vortex density for different r, in comparison with those in the dynamic 
relaxation starting from ordered and disordered states. The vortex density is defined as 

v(t,r) = (\v p \),v p = E - *i(*)]/2ir, (17) 

p(r,ip) 

where $i and 8j denote the orientational angles of Si and Sj, (0j — 9j) are valued within 
the interval [— 7r, 7r], the sum is over the four links of the clockwise plaquette at 
site (r,<f), and (■ ■ •) represents both the statistical average and the average in the ip 
direction. Numerical results are shown in Fig. [3j Outside the core of the vortex, e.g., 
at r = 127.5, the vortex density is initially zero, then increases with time, and finally 
reaches the steady value in equilibrium. This dynamic process is relatively fast, the 
same as that in the dynamic relaxation with the ordered initial state. In the core of 
the vortex, e.g., at r = 0.707 or 1.58, the vortex density initially increases rapidly to 
a large value, which even exceeds that of the disordered initial state, then decreases 
slowly and relaxes to the equilibrium. This dynamic process is even slower than that 
with the disordered initial state. Therefore it is not surprising that a strong logarithmic 
correction to scaling emerges. 
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Here we should emphasize that the slow dynamic relaxation in the core of the vortex 
state is indeed induced by the topological structure. For example, we may simulate 
the dynamic relaxation for the initial state with a four-spin vortex on the ordered or 
disordered background. The pseudo-magnetization rapidly drops to zero. 

To clarify the scaling form of the Binder cumulant for the dynamic relaxation of a 
vortex state, we plot U(t, r) as a function of t for different radius r in Fig. 11(a). Similar 
to Eq. (flE]) . we expect the scaling function U(s) in Eq. ffTO]) to be 

tt( \ r const s ~~* 00 CIQ\ 

U(s) ~ \ . . 18 

v ' 1 s~ do s^O K J 

In the core of the vortex state, e.g., r = 0.707, U(t,r) looks like exhibiting a power- 
law behavior, and the slope of the curve is 0.94(2). Assuming ~ t 1 ^, one derives 
U(t) ~ t&- x+ *M*. From (d-l + d )/z = 0.94(2), one calculates d = 0.88(4). With the 
logarithmic correction in Eq. ( fl2l) . however, the fitting yields do = 1.04(6), close to 1. To 
further confirm this result, we may fully examine the scaling form in Eqs. ( flOl) and ( TT8l) . 
In Fig. BJb), U(t, r) is plotted as a function of s = r/£(t). Data collapse is observed for 
different time t. In the core of the vortex, the logarithmic correction to scaling has been 
taken into account. Obviously, U(t,r) exhibits a power-law behavior in the core of the 
vortex. From the slope of the curves, one measures d—l + do = 1.97(4), then calculates 
d = 0.97(4), and it supports that d is close to 1. For comparison, the exponent d Q 
for the domain interface is close to 2 |22j. The reason is that the domain interface is a 
one-dimensional object, while the core of the vortex is effectively zero- dimensional. It is 
interesting that the pseudo-magnetization of the zero-dimensional vortex core exhibits 
a similar dynamic behavior as that of the one-dimensional domain interface, but its 
fluctuation is different. In fact, the fluctuation of the pseudo- magnetization at the core 
of the vortex state does not show r-dependence. Therefore, we have the simple relation 
d Q = rjo - (d - 1). 

Outside the core of the vortex, e.g., r = 100, a power-law behavior is observed in 
Fig. Eta), and the slope of the curve is 0.504(11). Assuming ~ one derives 
U(t,r) ~ t^" 1 )/ 2 . From (d - l)/z = 0.504(11), we obtain z = 1.98(4), consistent with 
the theoretical value z — 2. At the same time, we may also examine the behavior 
U(t,r) ~ in the large-s regime in Fig. Hfb). From the slope of the curves, one 

measures d — 1 = 0.996(9). Since the maximum radius of the vortex line is r = L/2 in 
the square lattice, the dynamic behavior of the Binder cumulant becomes unreasonable 
for r > L/2. Due to this boundary effect, the data collapse of the Binder cumulant 
in Fig. H(b) is less good for large s, compared with that for small s and for the 
magnetization shown in Fig. [2(b). 



3.2. Vortex state with quenched disorder 

To study the dynamic effect of quenched disorder, and to further demonstrate the 
importance of the topological structure of the vortex state, we randomly fix the 
orientations of the four spins at the center of the vortex, and repeat the simulation 
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of the dynamic relaxation of the vortex state and analysis of the pseudo-magnetization 
according to Eqs. (Q and ([TBI) . The scaling function M(r /£(£)) is plotted in Fig. E^b), 
and data collapse is also observed. However, the vortex exponent is modified to 
?7o/2 = 0.80(2), significantly different from r] /2 = 0.99(2) for the case without disorder. 
Such a scenario is similar to the dynamic relaxation around a disordered surface [2T1 13T] . 
But here disorder exists only at four lattice sites in the core of the vortex state. 

Topologically, the vortex state may be compared with a domain wall or a surface. 
It is interesting that the dynamic relaxation of the vortex state is similar to that of 
a domain wall rather than a surface, i.e., i] /2 = 1.0. However, a limited amount of 
quenched disorder around the center of the vortex could modify the dynamic universality 
class, leading to 770/2 = 0.80(2). This again shows that the topological structure of the 
vortex state is crucial. 

3.3. Spin-wave state 

For the spin-wave initial state, the time evolution of the pseudo-magnetization is 
displayed in Fig. Eta). Since the spin- wave state is locally ordered, the dynamic 
relaxation of the pseudo-magnetization is expected to be similar to that with the ordered 
initial state, M(t) ~ t~ r i/ 2z ( without an logarithmic correction to scaling. In fact, if 
one computes the vortex density v(t), it completely overlaps with that of the dynamic 
relaxation starting from an ordered state as shown in Fig. [3j In Fig. [5ta), the exponent 
7]/2z = 0.0593(5) is estimated from the slope of the curve of L x = 512 and 1024. A 
power-law correction to scaling refines the measurement to n/2z = 0.0586(3), with 6=1 
and c = —11.1. 

For the dynamic relaxation with both the vortex and spin-wave initial states, the 
standard magnetization is zero from the beginning of the time evolution. Only the local 
pseudo-magnetization could characterize the dynamic systems. 

For the spin-wave initial state, the time evolution of the Binder cumulant is 
displayed in Fig. E^b). The dynamic relaxation of the Binder cumulant is similar to 
that with an ordered initial state. Assuming £(£) ~ t 1 ^, one deduces U(t) ~ t*^ 1 )/ 2 . 
In Fig. \Mb), the slope of the curves of L x = 512 and 1024 is estimated to be 0.520(3). 
Obviously, there exists a correction to scaling in the growth law of £(t), and the fitting 
with Eq. ( TTTj) yields an exponent {d—l)/z — 0.507(8), consistent with the theoretical 
value z = 2. Finally, the finite-size dependence U(t) ~ can be also verified with 

the data in Fig. [5](b). 

4. Theoretical calculation 

Based on the linearized long- wave approximation, we may theoretically derive the scaling 
forms for the dynamic relaxation with the spin-wave and vortex initial state. Here we 
just report the main results, and leave the detailed calculations in the Appendix 1. For 
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the spin-wave initial state, the local pseudo-magnetization is independent of R, 



M p {t,R)(xr^ 2z . 



(19) 



This power-law behavior is in agreement with the scaling form in Eq. (TUl) . For the 
vortex initial state, the scaling form of the pseudo-magnetization is 



where r' = r/y4a 2 t/T oc r/t 1 /" = s. This result is principally in agreement with the 

scaling forms in Eqs. (Q and (fT6|) . although the logarithmic correction is still absent due 
to the linearized long- wave approximation in Eq. (121 p . For small r', M p (t, 

r >) ~ t-^ 2z s, 

therefore, rj /2 = 1. 
5. Conclusion 

With Monte Carlo methods, we have simulated the dynamic relaxation starting from 
both a vortex state and a spin-wave state at the KT phase transition of the 2D XY 
model. The dynamic scaling behavior of the local pseudo-magnetization and Binder 
cumulant is carefully analyzed. Although the core of the vortex state is effectively 
zero- dimensional, it exhibits a universal scaling behavior, controlled by both the bulk 
exponent 77 and the vortex exponent 7] . Taking into account the logarithmic correction, 
the vortex exponent r/o/2 = 0.99(2) is extracted. A limited amount of quenched disorder 
in the core of the vortex state alters the dynamic universality class, and the vortex 
exponent is estimated to be r] /2 = 0.80(2). In Table [U all the measurements of the 
critical exponents are summarized, in comparison with those in the literatures. Finally, 
theoretical calculations based on the long-wave approximation are presented, and the 
scaling forms for the dynamic relaxation with the spin-wave and vortex initial states are 
derived. 

Our numerical and theoretical results show that the local pseudo-magnetization is 
a useful concept in characterizing the relaxation dynamics starting from a metastable 
state or a low-energy excited state. Further applications of the methodology to complex 
dynamic systems such as spin glasses are interesting. 
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Appendix 1. Theoretical calculation 

In this appendix, we theoretically derive the scaling forms for the dynamic relaxation 
with the spin-wave initial state and the vortex initial state, based on the linearized 
long- wave approximation. In the Hamiltonian in Eq. ([Q), the spin Si can be represented 
by an angle, i.e., Si = (cos 9i, sin 9j). With the linearized long-wave approximation, the 
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Hamiltonian can be reduced, 

H = - £ cos(0, - S ) ~ | E - e if + H o 

<ij> * <ij> 

= ^E(W-^ + «)) 2 + ^o, (21) 

where R is the position vector of a lattice site and a is the spacing vector between 
the site and its nearest neighbors. The constant H will be ignored in the following 
calculations. After the fourier transformation, the Hamiltonian is rewritten as 

H= 1 -Y,J(k)\e(k)\ 2 . (22) 

k 

Here J(k) can be simplified with the long-wave approximation, 

J(jfe) = 1 E I 1 - e^f ~ a 2 /c 2 . (23) 

a 

The relaxation dynamics of model A for the XY model can be described by the 
Langevin equation, 

dB(t,k) 1 dH . r . a 2 k\ . - 

V = - 7^ — + eft, fe) = 7z-0{t, k) + e(t, k), 24 

alt Td6{t,k) K J T 1 ; ' 1 ; 

0(t,jfe) = JL / dRe~ ikli 9(t,R). (25) 

27T J 

Here e(t, fc) is a Gaussian white noise with (e(i, k)e(t' , k')) = 25(k + k')S(t — if) , and (■ • •) 
represents the statistical average. Since the equation is linear, it can be exactly solved, 

6(t, k) = f e-'^-WeW, k)dt' + 9(0, k)e- a2k2t ' T . (26) 
Jo 

One may easily calculate (9{t,k)) = 9(0, k)e- a2kH l T and 

(\9(t,k)\ 2 ) = (l- e-^/r) + \6(0,k)\ 2 e- 2a2k2t ' T . (27) 
or A; v 7 

In Eq. ( J26l) . 9(t,k) consists of two parts, the bulk one F(t,k) and the initial- 
condition dependent one G(t,k), 

F(t, k) = f e- a2k2(t - t,)/T e(t', k)dt' , G(t, k) = 9(0, k)e~ a2k2t/T . (28) 



From the inverse fourier transformation, one calculates F(t,R) and G(t,R), and then 
obtains the complex magnetization, 

M x (t, R) + iM y (t, R) = (e id ^) = (e lF ^) e iG ^ R \ (29) 

Using the cumulant expansion and the result (F(t,R)) = 0, one deduces 
M x (t, R) = e -(F 2 (t,Rm cos (G(t, R)), 

M y (t, R) = e -(F 2 (tM)/2 sin (G(t, R)). (30) 
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One may calculate 

(F 2 (t,R)) = (j^f J dk J dk' (F(t,k)F(t,k'))e< n+n ')- A 

= L\ nt + c, (31) 

where C is an integral constant. Obviously the bulk part of 6(t, R) is independent of R. 
Now it is important to calculate the initial-condition dependent part, 

G(t, R) = ^~ f dke^ n G(t, k). (32) 

Z7T J 

For the ordered initial state, i.e., 9(0, R) = (ft, one calculates G(t,R) = (ft, and 

M x (t, R) oc t~ T/87T cos 0, M y (t, R) oc t~ T/87T sin 0. (33) 

Under the long-wave approximation, n = T/2-k and z = 2 |46j. Thus the power-law 
behavior of the magnetization, i.e., M oc t~ v ^ 2z , agrees with that from the scaling 
theory and renormalization-group calculations [UE]. 

For the spin-wave initial state, i.e., 8(0, R) = 2<f> cos(k x) with R = (x,y), 

G(t,x) = — ( dke iin G(t, k) = 2<j)cos(k x)e- a2k ° t/T . (34) 
2n J 

Since we only consider the case of small k , e.g., k ~ 1/L, the time dependence in 
Eq. (I3"4"|) can be ignored. Hence G(t,x) is approximately equal to 9(0, R). Then 

M x (t,R) oc r v/2z cos(2(pcos(k x)), M y (t,R) oc r ,r,/2z sin(20cos(A;ox)).(35) 

Obviously the standard magnetization is zero. But the local pseudo-magnetization is 
non-zero. In fact the local pseudo-magnetization is independent of x, 

M p (t, R) = M(t, R) ■ M(0, R) oc r v/2z . (36) 

This power-law behavior is in agreement with the scaling form in Eq. (fT4l . 

How is the dynamic relaxation of a vortex state? In this case, 0(0, R) = <p(R) + ir/2, 
and ip(R) is the polar angle of R. According to Eqs. (1551) and (!32j) . one can calculate, 

G(t, r) = — ( dR9(0, R) e ~(.R~r) 2 /i^t/T) _ (37) 
Airt J 



Let us define R' = (R - r)/ ^Aa 2 t/T,r' = r/^AaH/T, then 

G(t, ?) = - J dR'R'e- R ' 2 J dip(R')9(R' + ?, 0). (38) 

Due to the rotational symmetry, G(t,r') — 9(0, r') is independent of f(r'). To simplify 
the calculations, we set r' = (r',0), i.e., (p(r') = 0, thus 

G(t,P) - 9(0,?) = - J dR'R'e~ R ' 2 J d<p{R')<p{R' + ?) 

= - dR'R'e~ / dif(R')Q(ip(R'),r'/R'), (39) 

7T JO J-tt 
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and Q(ip(R'),r' /R') is introduced for convenience, 

r' 

Q(v(R'),r'/R') = arctan[cos^( J R') + — , sm<p(R')}. (40) 

Here arctan(x, y) is defined in the interval [— 7r,7r] by tan[arctan(x, y)] = y/x. 

In fact, ip(R') is the polar angle of R', and Q((p(R'), r'/R') is that of R' + r'. r' 
is the fixed site we are looking at, and R' is the site we should integrate. For R' < r', 
\Q(ip{R'),r'/R')\ < tt/2, and B{(p(R'),r' / 'R') is a single-valued odd function of ip(R'). 
The integration over tp(R') in Eq. (1391) is zero. For R' > r', Q((p(R'),r' / R') is still 
odd, but not single-valued at <p(R') = ±tt. When <p(R') ±ir, Q(<p(R'),r' /R') -»■ ±vr. 
However, = ±7r correspond to a same spatial point. Since Q(ip(R'), r'/R') is 

multiple-valued at ip(R') = ±vr, the integration over ip(R') in Eq. (1391) can be non-zero. 

From the above analysis, the integration in Eq. (1391) can be rewritten as 

G{t,r')-e{Q,r') = - dR'R'e / d^(R')e(^(R'),r'/R'). (41) 

Let us focus on the integration over f(R')- Since Q(if(R'),r'/R') is an odd function, 
non-zero contribution only comes from a small interval around <p(R') = ±7r, and 

f ^[©(^ + e) - S(<p - e)] = 2 f ^9(^), (42) 

where e is a small angle to be determined. Substituting Eq. fflSl) into Eq. flHJ, 

G{t, r') - 6(0, r') = — / dR'R'e^ / <fy>[e(y> + e) - 9(^ - e)] 

27T Jr' 

= j_ rdR'R'e- R ' 2 r de4 

= ± A / dR'e- R . (43) 

Here 2e = A/R', and / dO = ±27r dependent on the existence of a vortex or anti- vortex. 

The pseudo-magnetization can be calculated by M p (t,r') = t"^ 22 cos(G(t, P) - 
0(0, r')). When r' — > oo, it approaches the bulk behavior M p (t,r') = t _r? / 2z . When 
r' 0, M p {t,r') = t^/ 22 cos(A-y/7r/2). Due to the symmetry, or from our numerical 
simulations, M p (t, 0) should be zero. Hence the parameter A = is determined. The 
finial form of the pseudo-magnetization is then derived 

MJt, r') = t^ /2z sin(v^ f dR'e~ R ' 2 ), (44) 

Jo 

where r 1 = r / ^Aa 2 t/T oc r jt 1//z = s. 

The linearized long-wave approximation in Eq. (F2"T|) is good for the dynamic 
relaxation of the spin-wave state, for the dynamic effect of the vortex-pair creation 
and annihilation is weak. For the dynamic relaxation of the vortex state, however, the 
linearized long-wave approximation could not catch the logarithmic correction in the 
core of the vortex. Additionally, a phenomenological technique has been used in the 
derivation of Eqs. ( 143]) and ( 1441) to obtain the scaling form in the core of the vortex. 
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STD [6] 


DW [22] 


Spin-wave 


Vortex 


Disorder 


V 


0.234(2) 


0.234(2) 


0.234(2) 


0.234(2) 




V2 




1.00(1) 




0.99(2) 


0.80(2) 


d 




1.99(1) 




0.97(4) 


0.78(3) 


z 


2.01(1) 


1.99(4) 


1.97(4) 


1.98(4) 





Table 1. Critical exponents estimated from the relaxation dynamics of the vortex state 
and spin-wave state, in comparison with those from the domain-wall (DW) dynamics 
and short-time dynamics (STD). In the last column, quenched disorder is located at 
the vortex center. 
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Figure 1. Dynamic relaxation of a vortex state and a spin- wave state is displayed 
respectively in (a) and (b), for the 2D XY model at the temperature T=0.89, slightly 
below Tkt- The arrowhead denotes the orientation of the spin. 
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Figure 2. (a) The time evolution of the pseudo-magnetization of the 2D XY model 
starting from the vortex state is displayed for different radius r on a double-log scale. 
Simulations are performed with L = 256, and solid dots are the data obtained with 
L = 512. Dashed lines represent power-law fits, and pluses indicate the fit with a 
logarithmic correction to scaling, (b) The scaling function M(r/£(t)) = M(i, r)£(i) T? / 2 
is plotted on a double-log scale. In the core of the vortex, the logarithmic correction 
to scaling is taken into account, and c = 5.45 is from the fitting in (a). Data collapse 
is observed for different t. The dashed line shows a power-law fit. 




Figure 3. The time evolution of the vortex density in the dynamic relaxation of a 
vortex state, in comparison with that in the dynamic relaxation starting from ordered 
and disordered states. 
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Figure 4. (a) The time evolution of the Binder cumulant of the 2D XY model 
starting from the vortex state is displayed for different radius r on a double-log scale. 
Simulations are performed with L — 256, and solid dots are the data obtained with 
L = 512. Dashed lines represent power-law fits, and pluses indicate the fit with a 
logarithmic correction to scaling, (b) Data collapse of the Binder cumulant is displayed 
for different t. In the core of the vortex, the logarithmic correction to scaling is taken 
into account, and c = 5.45. Dashed lines show power-law fits. 
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Figure 5. In (a) and (b), the pseudo-magnetization and Binder cumulant of the 2D 
XY model starting from the spin- wave state are plotted on a double- log scale. Dashed 
lines show power-law fits, while crosses indicate the fitting with a power-law correction 
to scaling. 



